Abstract. The semiclassical limit of the focusing Nonlinear (cubic) Schrodinger Equa- 
tion (NLS) corresponds to the singularly perturbed Zakharov Shabat (ZS) system that defines 
the direct and inverse scattering transforms (1ST). In this paper, we derive explicit expres- 
sions for the leading order terms of these transforms, which are called semiclassical limits 
of the direct and inverse scattering transforms. Thus, we establish an explicit connection 
between the decaying initial data of the form q(x, 0) = A{x)e lS ^ and the leading order term 
of its scattering data. This connection is expressed in terms of an integral transform that can 
be viewed as a complexified version of an Abel type transform. Our technique is not based 
on the WKB analysys of the ZS system, but on the inversion of the modulation equations 
that solve the inverse scattering problem in the leading order. The results are illustrated by 
a number of examples. 
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1 Introduction 

Direct and inverse scattering transforms play the key role in the solution of integrable sys- 
tems. Important results about the scattering transform with decaying initial data (potential) 
can be found in [26J. Our main interest lies in the semiclassical limit of the scattering trans- 
form for Zakharov - Shabat (ZS) system 

i —w ('- '< nr. it: 



dx \q -z 

where z is a spectral parameter, e is a small positive parameter and W is a 2 by 2 matrix- 
function. ZS system (CQ) is the first (spatial) equation of the Lax pair for the focusing 
Nonlinear Shrodinger equation (NLS) 

ied t q + -e 2 d 2 x q + \q\ 2 q = , (2) 

where iel and t > are space-time variables. 

The scattering data, corresponding to the initial data q(x, 0, e) consists of the reflection 
coefficient r init (z,e), as well as of the points of discrete spectrum, if any, together with 
their norming constants. Since the time evolution of the scattering data is simple and very 
well known ([25]), the evolution of a given potential can be obtained through the inverse 
scattering of the evolving scattering data. The inverse scattering problem for the NLS (j2J) 
at the point x, t can be cast as a matrix Riemann-Hilbert Problem (RHP) in the spectral 
z-plane, which is stated as: find a 2 x 2 matrix- valued function m(z) = m(z;x,t,e), which 
depends on the asymptotic parameter e and the external parameters x, t, such that: i) m(z) 
is analytic in C\r, where the contour r = R with the natural orientation; ii) 

fl + rf f\ T/ fQ . 

m + = m_ I I = m_ V {6) 
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on the contour T, where r(z,e) = r init (z, e) exp\^:(2z 2 t + zx)\ and m±(z) = lims_> m(z ± i5) 
with 5 > and z G R; iii) lim 2 _^ 00 m(^) = J, where J denotes the identity matrix. In the 
presence of solitons the contour T contains additional small circles around the eigenvalues 
with the corresponding jump-matrices (see, for example, [22J or [13]). 

In a more general setting, corresponding to AKNS systems, f(z) in the the jump matrix 
V should be replaced by p(z), which represents another piece of the scattering data that is 
independent of r(z); i.e., there is no functional dependence between r(z) and r*(z) on R; 
however, this case will not be considered in the present paper. It is well-known (see, for 
example [27]) that the RHP (J3]) has a unique solution m{z) that has asymptotics m{z) = I + 
— + 0(z~ 2 ) as z — > oo, and that the solution to the NLS (j2J) is given by q(x, t, e) = — 2(m 1 ) 12 , 
where (mi)x2 denotes the (1,2) entry of matrix m\. In the case when r init (z,e) has analytic 
continuation into the upper halfplane, the RHP for m(z) can be simplified by factorizing the 
jump matrix 



and "splitting" jump condition ([3]) into two jumps: one with triangular jump matrix V + 
along some contour T + in the upper halfplane C + (R is included in C + ) and the other with 
triangular jump matrix V_ along some contour T_ in the lower halfplane C~. Contours T-i- 
are deformations of R. Due to the Schwarz symmetry of ZS problem, contours T± can be 
choosen to be symmetrical to each other with respect to the real axis, and we can restrict 
our attention to only one jump condition, say, on the contour T + C C + . 

A contour T + G C + , which is a smooth deformation of R, together with a function f(z), 
which is analytic (or even Holder continuous) along T + with Qf(z) < —S for all suffisiently 
large z G r+, 5 > 0, define a solution q(x,t,e) of the NLS (j2J) in the following way: if rh(z) 
is the normed at z = oo solution of the matrix RHP with the jump matrix 

V + = (]. !jj , where r = e~^^ z) (5) 

on the contour r + , and the corresponding symmetrical jump V- = V£ (see 01])) on the 
symmetrical contour T_, than 

q(x,t, e) = — 2(mx)i2, where rh(z) — I -\ V 0(z~ 2 ) (6) 

as z —>■ oo. This construction holds even if f(z) depends on e. 

In the semiclassical limit problem (j5J), we consider initial data (potential) of the form 

q(x,0,e) = A{x)e lS{x)/£ , (7) 

where the amplitude A(x) is decaying at ±oo, and derivative of the phase S'(x) has the 
limiting behavior 

lim S'(x) = n± (8) 

x— >±oo 

with some finite m_ < In order to calculate the leading order of the solution 

q(x, t, e), t > 0, one needs to find the leading order of the solution m to the RHP ([3]) as 

£ -> 0. 
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The nonlinear steepest descent method ([9], [TO] ) , together with the ^-function mechanism 
([5]), is, perhaps, the most powerful tool of the asymptotic analysis of the RHP The 
key part of this method, in the setting of our problem (genus zero region), is the following 
scalar RHP for the unknown function g(z) = g(z; x,t), that: i) is analytic (in z) in C\ 7 m 
(including analyticity at oo); ii) satisfies the jump condition 

9+ + g- = fo - xz - 2tz 2 on 7 m , (9) 

for x G M and t > 0, and; iii) has the endpoint behavior 

3 

g(z) = 0(z — a)2 + analytic function in a vicinity of a. (10) 

Here: 7 m is a Schwarz-sym This observation remains true ifmetrical contour (called the 
main arc) with the endpoints a, a, oriented from a to a and intersecting K only at fi+; 
g± are the values of g on the positive (left) and negative (right) sides of j m , and; function 
fo = fo( z ), representing the scattering data, is Schwarz-symmetrical and Holder-continuous 
on 7 m . Taking into the account Schwarz symmetry, it is clear that behavior of g(z) at both 
endpoints a and a should be the same. Assuming f and 7 m are known, solution g to the 
RHP (Q without the endpoint condition ffTUl) can be obtained by Plemelj formula 

9(z) = m f 7F m),c , (id 

2m Jim (C - Z)R{0 + 

where 

f(z) = f(z; x, t) = f (z) -xz- 2tz 2 (12) 



and R(z) = a/ (z — a)(z — a). The branchcut of R coinsides with 7 m and the branch of R 
we use is defined by 

lim = -1 . (13) 

z— >oo z 

If fo{z) is analytic in some region S that contains 7 m \ {/i+}, the formula for g(z) can be 
rewritten as 

R(z) f /(C) , t hA s 

where 7 m C S is a negatively oriented loop around 7 m (which is "pinched" to 7 m in 
where / is not analytic) that does not contain z. 

It is well known (see, for example, [12J) that additional smoothness at the endpoints put 
some constrains on the location of these endpoints. To state these constrains, we introduce 
function h = 2g — f. According to flU 



h(z) = m f /(o ic , (15) 

where z is inside the loop j m . The endpoint condition (TTU|) can now be written as 

3 

h(z) — 0(z — a) a as 2 — >• a, (16) 
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or, equivalent ly, 

/(C) 



..(c-awo;*- - (17) 



the latter equation known as a modulation equation. Substituting (jTzJ) into (IT7I) yields 

x + 2(3?a + a ) t=7 L [ /o(C > ^ • (18) 
2 ™ Am (C - ")^(C)+ 

It is clear that for a given fo(z), (TT5|) defines a as a function a = a(x, £). The significance of 
a(x,t) is that it represents the leading e-order term 

q (x,t,e) = A(x,t)ei s ^ (19) 
of the solution q(x,t,e) to the Cauchy problem (E])-© through 

a(x, t) = a(x, t) + ib(x, t) = —-S x (x, t) + iA(x, t) . (20) 

By Schwarz symmetry, modulation equation (TTTI) holds if a is replaced by a. Adding and 
subtracting these two equations and using integration by parts, we obtain system of moment 
conditions 

||) rfc = , /K_£fK) dc = , (2!) 

<7m R \Q+ J Ira H \0 + 

where a = a + ib, which is equivalent to (fITj) . If fo is defined only on the contour 7 m 
(nonanalytic case), the loop integrals in (1211) should be replaced by the integrals over the 
contour j m . Substituting (fl"2l into (I2TT) . we obtain 

±fML dC = x + Ua, -L / K_^K) rfC = -2* . (22) 
2vr« ^ m i?(C)+ 2?rz i?(C)+ 

For a fixed t > 0, the second moment equation (j22j) defines a curve S in the spectral plane, 
whereas the first moment equation f!22p determines a parametrization of £ by x G R. 

Solving system fT22|) for a given /o( 2 ), i-e., finding a(x, t) that satisfies ff22l for all a; G R 
and all t G [0, to] with some to > 0, is the central part of the inverse scattering procedure 
for the leading order solution of the Cauchy problem d2])- (EJ) . Considerable progress has 
been achieved in solving this problem, see [22], [13], [21]. However, calculating fo, which 
represents the leading order term of the spectral data corresponding to (J7j), continue to pose 
a considerable challenge. In particular, the spectral data considered in [22] and [13J was 
calculated explicitly, because ZS system ([1]) for the corresponding initial data was reduced 
to the hypergeometric equation. In the case of general analytic initial data ([7]), the WKB 
analysis of singularly perturbed ZS systems in the complex x-plane seems to be the 
most natural approach for the direct scattering. There are, however, considerable difficulties 
associated with this approach even for relatively simple initial data, such as, for example, the 
need to keep track of a large number of turning points, singularities, Stokes lines that connect 
them, etc., (see, for examle, [H]). This is why, in our opinion, the results about the direct 
scattering are quite limited: one can mention numerical simulations of the discrete spectrum 
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of ([T]) in [2J , followed by formal WKB calculation in [16] of the F-shaped spectral curve from 
[2J, and rigorous WKB construction of discrete spectrum for certain special potentials ([7J) in 



The main goal of the present paper is to derive an explicit formula for fo(z) that will be 
valid for a rather broad class of initial data (jTJ) directly from the moment conditions (1211) 
(Section [2]). We proceed with studying properties of the transforms that connect fo(z) with 
the scattering data (Sections El - ED and, at the end, consider a number of examples that 
include already studied potentials (j7J), as well as some new cases (Section [7j). 

For the rest of the paper, unless specified otherwise, we assume t = 0. Then the first 
equation of (1221) can be considered as a transformation 

x(a) = — ! . f ^ (23) 
m 7 7m ^(C-a)(C-a) + 

of a given /q(C) into x(a), which has the meaning of the inverse function to a(x) = 
a(x, 0), x G R. Provided that a(x), determined by initial data ([Zj) through ( l20l) . is in- 
vertible for all x(a) is a real valued function defined on the curve £ that is the graph 

of a(x), i6l. The main result of this paper is the formula 



f Q (z) = / z - //+ + -^/(z - - w) a/(ii)du + (z - fJi+)x(z) + / (// + ), (24) 



where zeE and the integral is taken along S, which is the inversion of transformation (123|) . 
Here fo(fi+) is a free real parameter and the branch of the radical satisfies normalization (TTB"]) . 
Transformations fT23|) and (1241) resemble the pair of ^46eZ trans formants for axially symmetric 
functions, which is convenient to write in the form (pQ) 



f°° 1 /"°° r#- 

M(0 = -2/ N'( V )^fif^ed V and N(rj) = — M'(£) j== , (25) 



r/ 2 



where ^, 77 6 R. Indeed, if a G «M and E C iR, the radicals in ( 1231) and ( 1241) become \/£ 2 — a 2 
and \/ z 2 — u 2 respectively, and transforms (1231) and (1241) become a pair of Abel transforms 
(with some extra terms in (1241) required for convergence). However, if a G R and EcK, 
( 123]) becomes finite Hilbert transform (note that, due to Schwarz symmetry, /q(C) has a jump 
2^3/0 (C) along the real axis) on [a, //+]. Transformations ( 1231) and (1241) will be referred to as 
Abel-Hilbert (AH) or complexified Abel transformations defined on a contour EcC. 

A given initial data ([7j) determines a(x) = a(x, 0) by (T20l) . where A(x, 0) = A(a;) and 
^(x, 0) = S'(x). We consider the following three objects, defined by initial data (jTj): 1) 
solution q(x,t,e) of the Cauchy problem (T2J), (J7|); 2) solution q(x,t,e) to (J2J) defined through 
©-(El), where / = fo{z) — xz — 2tz 2 with fo(z) being the AH transformation of x(a) given 
by ( 1241) . and T + = S (here we assume that a(x), x G R, is invertible with the inverse a; (a); 
3) function q (x,t,e) that is defined by ( 1T9|) . ( 12U1) . where a(x, t) is a smooth in x G R and 
£ > solution of the moment conditions (122|) (analyticity of the initial data on 1 G R is 
required to define q (x,t,e) with t > 0). Since q(x, 0, e) = 9o(x, 0, e) (transformation (1231) is 
inverse to ( 1231) ). we call qo(x,t,e) a leading order semiclassical solution to Cauchy problem 
©, (J71), or simply a semiclassical solution. The main achievement of this paper is formal 
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construction of the semiclassial solution q (x,t,e) with some positive t for a given analytic 
initial data (J7J). Does the semiclassical solution q (x,t,e) indeed represent a leading order 
behavior of q(x, t, e) for t > 0? Although we do not have the general answer to this question, 
the following facts and observations indicate that the answer should be affirmative at least 
for some substantial class of the analytic initial data ([7j). 

Requirements guaranteeing that q (x,t,e) is 0(e) close to q(x,t,e) on a compact subset 
D of the x,t plane that do not contain breaking points (see below), can be found in [22J, 
[23]. The first requirement is that 

w(z) = sign(/i + — z)^sfo(z) < for z < \i- and for z > //+; (26) 

more precisely, w(z) is separated from zero if z < //_ and if z > //+ except for neighborhoods 
of where it has behavior 0(z — n±) respectively. Let 7 C be a bounded smooth oriented 
contour in the upper halfplane, called complementary arc, that connects a and /i_ and does 
not intersect 7 m (except at the common endpoint a). The second requirement is that for 
any (x,t) G D there exist main and complementary arcs 7 m and j c , connecting a(x,t) with 
/jl± respectively, such that the signs of $$h(z), where h is defined by ([15]) . satisfy inequlities 
(sign distributions) 

^sh(z) < on both sides of j m and Qh(z) > on at least one side of7 c . (27) 

As shown in Section[5l violation of smoothness of qo(x, t, e) in the process of time evolution 
leads to the break of the anzatz ( fl9l) for the semiclassical solution. Since ( |26i) is independent 
of x, t, the break at some point (xb,tb) means that at least one of the inequalities (I2T1) is 
violated at some point(s) z b G 7 m U j c , Zb ^ a(xb,tb), or that condition (fTBT) becomes 

3 

h(z) — o(z — a) 5 as z — > a, (28) 

if Zb = a(xb, tb). If Zb is not a branchpoint of fo(z) (regular break), then the situation can be 
corrected by introducing additional main and complementary arcs in the RHP fl9|)- frro|) . see 
[22] for details. This corresponds to the change of genus from zero to a positive even number 
(the genus is even because of Schwarz symmetry) of some hyperelliptic Riemann surface 
1Z(x,t) that shadows the evolution of q(x,t,e), see Remark ll.ll below. Correspondingly, 
the semiclassical solution qo(x,t,e) can be expressed in terms of Riemann theta functions 
defined by 1Z(x,t), and 0(e) closeness between q (x,t,e) and q(x,t,e) extends to the region 
beyond the break. 0(e) accurate approximation of q(x,t,e) by qo(x,t,e) can be extended 
through further breaks (see [21]) provided that the breaks are regular. The case when z^ is a 
branchpoint of fo(z) (singular break) requires additional study. Approximation of q(x,t,e) 
by qo(x, t, e) allows us to call qo(x, 0, e) the semiclassical limit of the initial data for solution 
q(x,t,e) (which is determined by fo(z))- 

Does go approximate solution q of the Cauchy problem (T5]), (j7])? The answer to this 
question depends on how accurately solution q approximates solutions q. A closely related 
question is how accurately the scattering data e~~f°^ (see (EJ) of q approximates the scat- 
tering data r init (z,e) of q. The authors are not aware of any general results of this nature, 
however, for a special family of initial data where explicit form of r init (z) is available (see 
[22]), 

f (z) = lim -ie In r ini t(z,e) (29) 

eraO Z 
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at every z in the domain of analyticity of lnr init (z,e) (with properly located branchcuts of 
fo(z)). Similar result was obtained in [13j for pure soliton solutions of (j2J) {r init = 0), where 
e ~Tfo{z) provided a good approximation to the corresponding discrete scattering data. The 
authors expect that (j29p holds for a wide class of general analytic initial data. That is why 
fo(z), obtained by AH transformation ff24|) . is called the semiclassical limit of the scattering 
data that corresponds to ©. The authors also expect that, subject to certain requirements, 
the semiclassical solution q Q is the leading order approximation of q as e — > 0. Establishing 
this fact seems to be the last remaining essential step towards the complete solution of the 
semiclassical asymptotic problem for the focusing NLS. 

Remark 1.1. The RHP can be modified to include N, N £ N, contours (main arcs) where 
g undergoes a jump. These contours define a hyperelliptic Riemann surface 1Z = lZ(x, t) of 
the genus N — 1, associated with the semiclassical limit of the initial value proble (Cauchy 
problem) d2])- (JTJ) . Higher genus of TZ indicates that the corresponding potential can be 
represented as a modulated iV-phase wave expressed through the Riemann theta-functions, 
see, for example, [22] • The authors belive that the AH transformation f[2~4"l) can be generalized 
to reperesent the semiclassical limit of the direct scattering transform for the higher genus 
cases, however, this paper is restricted to study potentials of the form ((Tj) only, i.e., to genus 
zero potentials. 



2 Inversion formula for the AH transform 

In this section we prove that under the appropriate assumptions on a(x) the transformation 
inverts the transformation (1241) . Let us assume that: 



1. a(x) = —~S'(x) +iA(x) is a complex valued C 1 function on R, where A(x) is positive 
and S'(x) satisfies flS]); 

2. a(x) is locally and globally invertible, i.e., a'(x) ^ on R and the graph S of a(x) 
does not have points of self-intersection; 

3. the inverse function x(a), a£E, satisfies 

lim (« — fj, + )x(u) = 0, uGS, and (u — fi + )x'(u) £ L 1 (S lt() ) (30) 

for any u £ S, where £ uo denotes the arc of £ connecting u and 

Theorem 2.1. Ifa(x) satisfies conditions® then fo(z) in (|24|) and its derivative are well de- 
fined. Moreover, transformation (|23|) is inverse to transformation (1241 . i.e., the substitution 
of (|24~|) into (|23|) turns the latter one into the identity. 

Proof. According to our choice (fl3|) of the branch of the radical in (124")) . we have 

V(2 - u)(z — u) ~ —(z - fi + ) (31) 
as u — > provided z is separeted from //+. Then 



z - //+ + V(2 - «)(^ — u) — — j n — , (32) 



2(2 — /i + )(/i + — + 


ji + — u\ 2 




2 -//+ - yj{z ~ U)(Z - U) 


y/(z ~ U){z - U) 
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where, according to (!3T|) . the denominator approaches — 2(z — /i+) 2 as«-» Now, conver- 
gence of the integral in (I24p follows from (I32p and the second condition of (I30p . Differentiation 
of ([MD yields 



1 + 



2 — Km 



yj(z-1l)(z-u) 

where the integrand can be expressed as 

z — 



x'(u)du + x(z) 



y/{z- u){z - u) 



(33) 



(34) 



Z — fl + — \J (z — u){z — U) y/(z — U)(z — u) 

and the second 



Since \^su\ < \u — fi + \, convergence of the integral in (1331) follows from 
condition of (13"U1) . 

Substituting of fl33|) into ff23l) . which can be converted into 



xia) = — c d 



2, 

7T 



/SCO 



^(C-aJCC-a). 



(35) 



yields 



2 ^ Q 
x(a) = — Q 

7T 



x /(z-«)(z-«) , 

x (ujcra + 



*(C) 



C ^(C-aXC-a) 



V(C-«)(C-«) 



2 

( 

7T 



3 



-Km 



•^/(z-u)(z-u) 



X[U) 



y/{C ~ «)(C - a) " \/(tt - a) (it - a) 



Denote by 



1 + 



p(C) 



(36) 



(37) 



v/(C-«)(C-«) 

and by £ a negatively oriented loop that contains contour E Q U E Q . Then j^p(()d( = 
by Residue Theorem. If S denotes a negatively oriented loop around contour S u U £„, 
and S-t denote positively oriented loops around contour S a \ S u and its complex conjugate 
respectively, see Figure [fl then 



So 



KCR+ p(C)dC- 



(38) 



=, since the remaining term of p(() attains the same 



However, <L p(C)dC = <fc - , . 

values on both sides of the branchcut S u U S u . Taking into the account Schwarz symmetry 
of the intgrands, we obtain 



3 / P (C)dC = 3 / 



m+ V(C-a)(C-a) 
9 



(39) 



Figure 1: Contours £, Sq and E-j-. 



Returning to fl36|) and taking into the account fl39|) and the that dx(u) = x'(u)du e 
when m 6 S, we obtain 



3 



x'(u) / p(()d( 



du 



x'(u)Q 



p(()d( 



du 



x'(u)^s 



d( 



/-<+ 



y/(C-a)(C-a) 



du =Q ; / x'{u) 



d( 



v/(C-a)(C-«) 



. (40) 



Thus, 



x(a 



- — 

7T 



2 

7T 



xIm)^ 



d( 



d( 



+ 



/' t- y/(C ~ «)(C - a) y/(u-a)(u-a) 



du 



y/(C-a)(C-a) 



x(a) 



1 



§ y/{C-a){C-a) 



x(a) , (41) 



where we used the fact that, according to the limit in (1301) . the contribution from u = //+ is 
zero. □ 



3 Derivation of transformation (EH) for /o(^) 



So far, the observation that transformation ( 1231) resembles Abel transformation helped us 
to guess transformation (1241) . to which (123]) is inverse. In this section, we will show how 



10 



transformation (|24|) can be derived from the analysis of the RHP (Q. In particular, we 
show that solution g(z) for the RHP represented by integral (TT4")) in the complex z-plane 
(Plemelj formula), can also be represented by a dual integral in the complex x-pane (see 
fl62|) below). The "input" data for the integral representation in the z-plane is fo{z) and 
x, whereas the "input" data for the dual integral representation in the x-plane is a(x) and 
z. Function h(z) has similar integral representations. Then teh semiclassical limit of the 
spectral data fo(z) is given by f = 2g — h + xz. 

We start our derivation with the following observation. 

Proposition 3.1. Conditions fl9l)- ffl0l) imply that = 0. 

Proof. Applying 



dR(z) R(z) 



and 



1 i 



(C-*)(C-a) z-a 



da 2{z — a) 

to ffTTT) . we obtain 

<!</(-) R(z) f /(C) 



(42) 



_j i_ 

C-2 C-« 



-dC- (43) 



<9a 8yri(z - a) J^ m (( - a)R(Q. 
But, according to fTT7|) . the integral in (T431 is zero. The proof is completed. □ 



Under our convention t = 0, functions R,g,h,f depends on x and z. Since z will 
be considered as a parameter for the rest of the paper (unless specified otherwise), it is 
convenient for us henceforth to put the variable x in these functions in the first position. For 
example, the radical R introduced in (flTT) can be now rewritten as 



R{x, z) = ^{z-a{x)) 2 + b 2 {x) (44) 

with the same choice of the branchcut in the z plane as before. The choice of the branchcut 
for R(x, z) in the complex x-plane is discussed below. 

Since f(x,z) = fo(z) — xz, Proposition 13. II implies that total derivatives 

-^g(x,z) = ^-g(x,z), and -^-h(x, z) = -^-h(x, z) (45) 

ax ox ax ox 

coincide with the corresponding partial derivatives. Differentiating both sides of the RHP 
(J9]) in x and using (fTUj) . we obtain 

dh(x,z) dg{x,z) 1 

dx =~R(x,z), dx =--[z + R(x,z)] (46) 

for all x £ R and z G C + , where z is considered as a parameter and C denotes the upper and 
the lower halfplane respectively (we also use notation B = B Pi C =t for any set B). The fact 
that derivatives in fT46l) are independent of / opens the way to reconstruct h(x,z), g(x,z) 
and, thus, fo(z). 

To construct fo(z), we require that 

a(x) = a(x) + ib(x) = ~^S'(x) + iA(x) (47) 

satisfies the following conditions (A) : 
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1. 



a(x) and b(x) are real analytic on R and b(x) > 0; 



2. all but finitely many points of the parametric curve 

OS = {z G C : z = -S'(x)/2 ± iA(x), (48) 

are regular points, i.e., the tangent vector a'(x) / 0, i 6 1; moreover, dS does not 
have points of self-intersection (note that we used notation £ for d£ + in Sections [TJ 
ED; 

3. if T> denotes the common domain of analyticity of a(x), b(x), then a(T>) D C; here a(x) 
is defined by (l47j) and £ is an open domain in the upper z-halfplane that contains the 
union of the strip {z : < Qz < L}, where L > 0, with S + , where £ denotes the open 
region of C bounded by d£; 

4. there exist ^£8 and p > 1 such that 

lim x p [a(x) — fi±] = lim x p b(x) = 0, (49) 

5Rx— >±oo iRx— >±oo 

respectively for all a; such that a(x) G C; 

The curve 9£ connects fi + and /i_ in both C + and C _ . Condition A2 implies that there 
exists the inverse function x(z), where x(a(x)) = x Vx G M. Note that x(z) is analytic on 
at all the regular points. Any point x* G V, such that a'(x*) = 0, is called a point of 
ramification. The corresponding z is called a logarithmic (log) point. We can analytically 
continue x(z) from dS + to £ with the exception of the log points. To define x(z) uniquely 
in, £ we make branchcuts connecting every log point z* G £, Qz > with K. U oo in such 
a way that they do not intersect £ + and keep C connected. Let C denote C with the cuts. 
Then B = x(C) is the image of C by the conformal map x(z), see Fig. [2j 

According to (|46j) . we have 

h(x,z) = - R(y,z)dy (50) 

Jx(z) 

for all z G C and all i G B, where the contour of integration lies in B and does not cross 
the branchcut T(z) of R(x, z) in the complex x-plane (z is fixed), which is a path connecting 
x(z) with — oo (see Fig. [3]). 

For any fixed z G C, function /i(x, 2) is analytic in x for all finite x G £> except the 
branchpoints x = x(z) and x = x(z) of the radical R (the latter may or may not be in B. 
On the other hand, 

dh^ = _ r x z - a(y) dy + y/(2:-a(x(z)))(z-a(x(0))) 

^ ' A(^) y/(z - a(y))(z - a(y)) a'{x(z)) 

where 



<5(x) = a(x) = a(x) — i6(x) . (52) 
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Figure 2: Map x(z) maps dS + into the real x-axis, £ + (darker area) into B and the rest of 
C (lighter area) into B + . Possible cuts are not shown here.) 




x(z) 



Figure 3: The branchcut of R(x, z) is shown by the dashed lines: a path connecting a(x) 
and x is fixed (left); a path connecting x and — oo, z is fixed (right). 

Since a(x(z)) = z in B, it is clear that the latter term is zero for all z except when a'{x(z)) = 
0, i.e., except when z = fi± or z is a log point. So, 

,) = _/" - -'- a ^> (53) 
^ y x (*) v(z - ~ 

if a' (#(2)) 7^ 0. That means that for any finite x £ B function h(x, z) is analytic in z G C 
except for the log points, the branchpoints z = a(x), z = a(x), and, possibly, points z = ji±. 

Does h(x, z) from (J50|) coincide with ft, given by (IT5|) ? If the corresponding a(x) and /o( 2; ) 
satisfy modulation equation (|17|) . then they may differ by an independent of x constant. 
However, this constant is identically zero since h(x(z),z) = for all z G C in both cases of 
(|50|) and (fl5]l . Thus, for any fixed x E B, 

3 

ft(x, ,2) = 0(z — a(x)) 5 as 3 — > a(x), (54) 

provided that a'(x) 7^ 0. Equation (j54j) also follows directly from f)53p . Finally, for all x G M, 
we extend ft(x, z) into the lower z-halfplane C~ by Schwarz reflection. Note that h(x, z) has 
a jump 2i^sh(x, z) for z G R. In the case z G C - , the values of h(x, z) can be analytically 
continued from iGlto the complex x-plane. 

To calculate g(x, z), we first note that, according to (|4"4|) . ffT31 and condition A3, 

R(x, z) = — (z — n±) + o (x~ p ) as — > ±00, x E T> (55) 
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for any fixed z G C. Then 



z — fi± + R(x, z) 



2(z - fi±)[a(x) - fi±] - [a(x) - fi±] 2 - b 2 (x) 



a[x) — fx± 



z — fj,± — R(x, z 
R 2 (x, n± 



+ o(x- 2p ) 



(56) 



as 3fcr — > ±00, x G T>, for any fixed z 7^ fi± respectively. According to (|46|) . we define g(x, z 
as 

1 



g(x,z) 



/ [z - n+ + R(y,z))dy + /i+x 



+ K(z) 



(57) 



for any finite x G B and z G C, where K(z) does not depend on x and the contour of 
integration is in B. Convergence of the integral in (jSTj) follows from fl56|) . 

Assuming that a(x) and /o(^) satisfy modulation equation ffT7|) . we want to determine 
so that g(x, z) defined by f|5T|) coincide with defined by ffTTj) . Note that 



lim 



#(x,z) + ~/i+x 



Rewriting ffTTl) as 



g(x,z) 



-R(x, z) 



JT(z) . 
/o(C) - *C 



4m L (C-^)i2+(a;,C) 



and taking limit of (159]) as x — > +00, x G R, x G £>, we obtain, according to 



lim 

B— »+00 



lim 

x — >~\-oo 



> - ^+)(/o(^+) -^+) 
47ri(z — /!+) 



*» t/(C-a(z))(C-a(x)) 



+ 2^+3 



;/o(a*h 



Comaring (158]) and (1601) . we obtain 



(58) 



(59) 



(60) 



(61) 



the latter follows from the requirement that g(x,z), x G R , is Schwarz-symmetrical in z. 
Thus, we obtain 



g{x,z) 



/X 
(z- /!+ + R(y, z)) dy + /i+x 
.-00 



(62) 



where /o(/ i +) is a f ree rea l parameter. 

We want to emphasize that (|62|) represents a new form of solution to the RHP (!9l- (jT0l) 
(with t = 0), written as an integral in the x-plane, see Theorem 13.21 below. 
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Note that for a finite fixed z G C + function g(x,z), defined by (I62I) . is analytic for all 
finite x G B except x = x(z) and x = x(z). According to 0551) . we have 

1 , z - a(x) a(x) - ^ + o{x _ p) 



R(x, z) Z — pL± 

as $tx — > ±00, x G V, for any fixed z 7^ ji± respectively. Thus, for any finite x G B, function 
g(x,z), defined by fl62|) . is analytic in all finite z G C except z = /x+ and branchpoints 
z = a(x), z = a(x). 

Now, using / = 2g — h, we obtain 

p+00 

f(x,z) = [z-fi + + R(y,z)]dy+(z-fi + )x(z)-xz + f (jj l+ ), (64) 

Jx(z) 

so that 

r+00 

fo(z) = / [«~A*+ + ^)]%+(2-A t +)a:(2) + /o(A t +)- (65) 

Jx(z) 

Since /i(x, z) and g(x, z) are analytic in z G C except the branchpoints z = a (a;), z = a(x), 
the log points and, possibly, points z — fi±, function fo(z) = 2g(x, z) — h(x, z) —xz is analytic 
over the same domain. But fo(z) does not depend on x, hence it is analytic everywhere in 
C except the log points and possibly, points z — with 



fo(z) 



x(z) 



' z - a(y) 

R(v,z)\ 



dy + x(z). (66) 



Theorem 3.2. For any a(x) satisfying conditions Al - A4, functions g{x,z) and fo(z), 
given by ( |62l) and ( |65i) respectively, satisfy the RHP (l9i)-(fT0l) aHieM swc/i i/iat a'(x) 7^ 0. 

Proof. The contour 7 m in the RHP ([9]) can be deformed within the domain of analyticity of 
fo (endpoints a(x),a(x) and the midpoint fi + remain fixed) without affecting the solution 
g(x,z). There are no more than a finitely many log points in any compact subset of L. 
According to the construction, set C is connected, i.e., without any loss of generality we can 
assume that 7 m lies within the domain of analyticity of fo (with the exception of the point 
/JL+). Let us fix some x G R, such that a!{x) 7^ 0, choose some z on the corresponding j m , 
and consider g+(z) + g~(z). The contours of integration for g± in fl62|) lie on opposite sides 
of the branchcut of R(x, z) in the complex x-plane (z is fixed), as shown on Fig. HI Then 



f+oa r x { z ) 

g+(x,z)+g-(x,z)= [z - //+ + R(y, z)]dy+ (z- n+)dy - fi+x + / (/i+) = f(x, z) . 

J x{z) Jx 

(67) 

To prove that g(x,z) is analytic at z = 00 for any x G B \ {x(oo)}, we fix some x and 
consider 

a? ( *' l) = -2/ w ( 1+ *tof)*' (68) 

where the contour of integration does not pass through x(oo). Then 

J2(y, z) = -(z - o(y)) + O f-^rl (69) 
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X 




Figure 4: The branchcut of R in B, connecting x(z) and — oo is shown by the dashed line, 
contours of integration for g± are shown by solid lines. 



as z — > oo uniformly on the contour of integration. Therefore, the integrand in fl68l) is of the 
order 0(b 2 (y)) uniformly on the contour of integration as z — ► oo, so that -§^g(x, z)\ z=OQ is 
well defined. Proof of fflQj) follows from (J54"|) and the fact that /o(^) is analytic at z = a (a;). 
Thus, requirements i) - iii) of the RHP (J9])- (JT0"]) are satisfied. □ 

Remark 3.3. Theorem 13.21 remains true if in the expressions for g(x, z) and fo{z) we replace 
+oo with — oo and /i + with 

Remark 3.4. For any z G S + , the contour of integration in fl64l) is the interval [x(z),oo) of 
the real axis. Then the change of variables u = a(y) converts (1641 into the transformation 
( l24l) . Note that this is not true for z G" £ + , since u ^ a(y) for complex y. 



4 Inequalities (ED-© for f 



In Sections |2] and [3] we constructed the semiclassical limit of the scattering data fo(z) for ZS 
system (pQ) with a given potential (JTj). Let us assume that potential © satisfied conditions 
A. In this section we discuss inequalities (126]) for /o as well as existence of the main and the 
complementary arcs satisfying (|27|) for all iGl. 

We start with the observation that for any x G K we have 

S/i(x,//±) = -9f / ^ (// ± - a(y))(fi± - a{y))dy = , (70) 

«/±oo 

since the contour of integration lies on the real line and the integrand is real valued. Con- 
vergence of the integral in ( 1701) follows from ( 1491) . 

Let us consider w(z),z G R. According to ( 1261) . ( 1501) and Schwarz symmetry of g(x,z), 
x G R, we have 



w(z) = sign(/i + - z)9 / a/(z - a(y))(z- a(y))dy, (71) 

where x G R can be choosen arbitrarily, i.e., the right hand side of (fTTjl does not depend on 
a particular choice of x G R. Therefore, we can choose contour of integration in ( J7TT) as the 
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Figure 5: Two solid curves in the x-plane show the boundaries of B. The image of the 
segment [x(z), x(z)] in z-plane is shown by the dashed line. 

vertical segment [x(z), x], where x = 3far(z), traversed in the negative direction, i.e., down 
(note that Qx(z) > 0). 

Lemma 4.1. If 

| arg(z -a(y))- arg(z - a(y))\ < vr (72) 
for every y G [x(z),8lx(z)], where z > //+, then w(z) < 0. If 

| arg(z — a(y)) — arg(z — a(y)) — 2n\ < n (73) 

for every y G [x(z),8lx(z)], where z < then w(z) < 0. 

Proof. Consider, for example, the case z > fi + . Since dy is negative purely imaginary, it is 
sufficient to show that 3? a/ (z — a(y))(z — a(y)) < 0. Taking into account our determination 
of the square root (jTBl . the latter condition is equivalent to | axg(z—a(y))+a,rg(z—a(y))\ < tt, 
which, together with (|52|) . implies (!72|) . see Fig. [51 Similar arguments prove the remaining 
case z < //_. □ 

Remark 4.2. Listed below are some sufficient conditions for (1721) and (1731 to be true for some 
z > /i + or some z < /i_ respectively. 

1. w(z) < for some z > fi + or some z < /i_ provided 

[^),^)]C8. (74) 

Indeed, this condition implies that and a(y) G for any y G [x(z),x], see Fig. [51 Since 
a(y) G C + \ £, the inequality (1721) is satisfied. Similar considerations are valid for (1731) . 
Condition (1741) is satisfied if, for example, the upper boundary <9,B + satisfies the vertical line 
test and if the complex conjugate dB + C B. 

2. According to lemma |4~1~1 w(z) < for some z > fi+ provided that for every y G [x, x(z)} 
the angle between z — a(y) and z — a(f7) is less than 7T. This happens, for example, if one 
can draw a line / through the point z so that the curve a(y), y G [x(z), x(z)] does not cross 
/ (lies in one the halfplanes produced by /). In particular, if / is a vertical line, we obtain 

®a(y) < z (75) 
respectively for all y G [x(z),x(z)]. Similar results hold for z < 
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Figure 6: Vectors z — a(y) and z — ac(y) for y £ [x(z), x] C M. 



We now consider inequality (127)) under additional assumption that a(x) is a monotonically 
increasing function. We choose 7+ to be the part of d£ + connecting a(x) and Take any 
z £ 7+. Then x(z) is real and x(z) < x. Let us choose any y £ z], it is easy to check, 

(see Fig. [6]) that 

< arg ((* - a{y)){z - ^))) < 2tt , (76) 

so that, taking into account (|5D)) and (1T31) . we obtain 5sh(x, z) > on 7 C . 

Let us now study the sign of 5sh(x, z) on the main arc 7+. Notice that h + + h- = 
on 7 m , so that, in general, Q7i has opposite signs on opposite sides of j m . Therefore, we 
have to have $sh = on j m . Let us first find sign^/i on the arc of d£ + between a(x) and 
fi + . Equivalently, we can think of deforming 7+ to the above mentioned arc of d£ + and 
finding signS/i on the negative (external) side of the arc. It is easy to see that if point 
z will pass over a(x) outside (of £), then the sum of the angles, shown at Fig. [6j i.e., 



arg y(z — a(y))(z — a(y))j , satisfies inequalities (176)) . Repeating the previous arguments, 

we readily obtain 5sh(x, z) < 0, z £ d£ + between a(x) and fi + on the outer (negative) side 
of the contour. Of course, on the opposite (positive) side of the branchcut, 5sh(x, z) > 0. 
That means, that inequalities on 7+ will be satisfied if there exists a branch of the level 
curve 5s h = (x £ R is fixed), connecting \i + and ot(x). This question is discussed below. 

We now want to calculate §j(#, /■*+) by taking limit z — ► /x+ along the negative (outer) 
side of 7+ (it is still asumed that 7 m deformed to coincide with arc of d£ + connecting a(x) 
and According to ()53l) . 



a/i. . z-a{y) 



dz 

Therefore, the fact that 



x,z)= / ; . (77) 



yj(z - a(y))(z - a(y)) 



< ~(y>i + ^2) < vr , (7* 



where y> = arg(z — a(y)), <p\ = arg(z — a(y)), y?2 = arg(z — a(y)), together with the 
determination of the proper branch of the radical R(y, z), imply S||(x, z) < in ( 1771 ). 
To prove (1751) . we first notice that the monotonicity of a(x) on K implies 

0<y?<|, -|<</?i<|, < v? 2 < |, ^2>^>ki|, (79) 
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<x(y) 



«(y) 




• z 


» 




/ H + 



cp 2 



a (y) 



Figure 7: Angles <p\ (negative), y?2 and cp for a(y), where x < y < x(z) 



see Fig. [71 so that the second inequality of ("78]) follows from (1751) . Notice that (3i > /3 2 , 
where Pi = | + /?2 = f — y?2- Then the remaining inequality (1751) becomes 



2tp>/3 1 -/3 2 . 



Let us inscribe the triangle a(y), z, a(y) into the circle. As shown on Fig. [HI cases (pi < 
and v?i > 0, both angles (5\ — /3 2 and 2y? rest on the arc z, O, z. Then (1501) follows from the 
fact that the vertex of angle /3i — is on the circle, whereas the vertex of angle 2<p is inside 
the circle. Thus, (1781) is proven. 



*1) 



The fact that lim z ^+ 5s^(x,z), if exists, is negative implies 



w'(/j,: 



lim S 



dh 

dz 



(x,z) > 0'. 



Choosing the other branch of R(y, z), we obtain 



lim ij— — (a;, z) < (J 
dz 



(82) 



if the limit exists, where /x+ is on the positive (inetrior) part of j m . 

Equations (1701) . (15*11) and (152]) show that 5sh(x, z) is negative for real z in a vicinity of 
Thus, there exists a zero level curve A of 5sh(x, z) emanating from z = //+ into the upper 
half-plane. Suppose that A C C and connects /i + with a(x). Then the signs of 5s h to the 
right and to the left of A have to be negative, i.e., the first inequality in (l2~7j) is satisfied. The 
results of this section can be summarized by the following statement. 

Statement 4.3. Let initial data ("TJ) be such that: assumptions (A) are satisfied and S'(x) 
is monotone on R; w(z) satisfies inequalities (1261) . If for a given x G R a zero level curve 
X, emanating from passes through a(x) and does not intersect branchcuts of f , then 
qo(x,0,e) is 0(e) approximation of q(x,0,e). 
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Figure 8: Triangle a(y), z, a(y) inscribed in a circle. The left circle shows the case <pi < 0, 
the right circle shows the case <p\ > 0. 

5 Some geometrical aspects of transition to a higher 
genus (breaking) 

In this section we consider t assuming nonnegtive values. Functions / and h are defined by 
( 1121) and ( fl5l) respectively. As it was established in [22], [21], transition from the genus zero 
to a higher genus (genus 2) occurs when an additional branch of zero level curve of Qh(x, t, z) 
(the branch that begins and ends at z — oo) intersects with the contour 7+ = 7+U7+ at some 
point zi, = Zf,(x,t), which is not a log point (the case of g + intersecting a log point requires 
further investigation). There are at least four zero level curves of ^sh(x, t, z) passing through 
the point z = Zb(x,t), which implies that h z (x,t,Zb) = 0. In the case Zb(x,t) 7^ a(x,t), 
the point 2& is called a double point. There are exactly four zero level curves of ^sh(x,t,z) 
passing through a double point Zb (degenerate cases of multiple level curves collision are 
not considered here, however, see [21]). In the remaining case Zb(x,t) = a(x,t), a triple 
point. There are five zero level curves of $sh(x,z,t) passing through a nondegenerate triple 
point Zb(x,t). The breaking point Xj,, t& on the x, t plane that corresponds to a triple point 
Zb is the starting point of the breaking curve that separates the genus zero and genus two 
regions for t > i&. Typically, the breaking curve forms a corner at Xb,h- A regular point of 
the breaking corve corresponds to a double point z&. The simple modulated wave qo(x,t,e) 
(semiclassical solution), given by (fl9i) . ( J20l) fails to approximate solution q(x, t, e) in the genus 
two region beyond t = t&. In this region q can be approximated by a two-phase modulated 
wave constructed through Riemann theta functions (see [22]). The onset of a two-phase 
wave (genus two) behavior of q corresponds to the appearance of a triple point a(xb,t b ) on 
the spectral plane. The main result of this section is that at the first breaking point x^, tb 
the semiclassical solution q (x,t,e) loses its smoothness. More precisely, if a(x, t) is a triple 
point for some x = Xb,t = t b then the derivative a x (xb,t b ) = 00. In other words, genus zero 
anzatz qo(x, t, e) experiences the first break at t — tb only if the amplitude of go or derivative 
of its phase (or both) develops an infinite slope at x = x&. The proof uses representation 
(|50|) of h(x, z) = h(x, tb, z), where tb is fixed, as an integral in the complex x-plane. 
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Lemma 5.1. If assumptions A for a(x) = a(x, t^) are satisfied then for any fixed xsKw 
have 

^-h(x,z)= y/z - a(x) + 0(z - a(x)) 3/2 (83) 

oz \/2a'(x) 

in a vicinity of the branchpoint z = a(x), provided a(x) is not a log point. 

Proof. Since a(x) is not a log point, we have 

dx 5 
x(a(x) + 5) = x + —5 + o(S) =x+ — — + o(5), 
dz a{x) 



where 5 G C is a small. Using (J53J), we calculate 



h z (x,z) 



R(x, z) 



. "^m z-a(y) 

— hm — r- / — . ==ay = 

4x) 5-> R(x,a(x) + 8)J x s/{z-a{y)){z-ct{y)) 

lim ; w v r = — , 84 

s^o a'(x)(2ib{x) + 5) 2a' (x) v ; 

which implies (1831) . □ 



As an immediate consequence of Lemma 15.11 we obtain 



h(x, z) 



v / 2tb(x) , — 

+ C(V-2 - a{x)) 



3a' (x) 



a(x)) 3/2 . (85) 



It is well known that zero genus anzatz for zero dispersion limit of the KdV equation 
breaks down when it develops infinite slope. The following corollary, which is another im- 
mediate consequence of Lemma 15.1} is an analog of this statement for the focusing NLS 
equation. 

Corollary 5.2. In the conditions of Lemma I5.il a(xh) = a(xb,tb) is a triple point only if 
a'(xb) = oo. 

Proof. a(xb) is a triple point only if the leading term of (|85|) is zero. Since b(x) > on IR 
(condition Al), we conclude that a'(xb) = oo. □ 

Remark 5.3. According to [5J, solution q(x,t,e) to Cauchy problem (T5|), ([7]) develops elliptic 
umbilic singularity at (xb,tb) if a(xb,tb) is a triple point. It is hypotesized there that in a 
vicinity of (xb, h) a solution q(x, t, e) can be approximated by the special tritronqee solution 
of the first Painleve equation PI (see [5] for further details). 

An obvious geometrical interpretation of Corollary (15.21) is the curve a(x) in R x C, where 
x G R and a(x) G C. For any x G R, this curve has a unique intersection with the plane 
perpendicular to R at x, which we interpret as a spectral plane C, see Fig. [91 This property 
is preserved under the NLS evolution a(x,t) = a(x,t) +ib(x,t), where a(x, 0) = a(x), as 
long as (x, t) is in the genus zero region. To get to a higher genus region, the curve a(x, t) 
must develop a "fold" (at least in the solitonless case), shown at Fig. [101 It is clear from 
the topological point of view that a(x,t) should become tangential to the spectral plane at 
some (xb,t b ) before the fold can develop, see Fig. [TTJ The point (xb,t b ) correspond to a 
triple point. Geometricaly, it is clear that a x (xo,to) = oo, which is exactly the statement of 
Corollary (E2J). 
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Figure 9: Curve a(x,t) in 1 x C, spectral plane C is orthogonal to real x-axis. Genus zero 
case. 





Figure 10: Curve a(x,t) in 1 x C, spectral plane C is orthogonal to real x-axis. The curve 
developed a "fold". Three point of intersection of a(x, t) and the spectral plane are indicated. 
For the values of x, corresponding to the "fold" , points (x, t) belong to the genus two region. 







oc(x,t) 






X 



Figure 11: Curve a(x,t) in 1 x C, spectral plane C is orthogonal to real x-axis. The point 
of tangency between a(x,t) and C corresponds to a triple point. 
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6 Symmetry 

In this section we establish connections between the symmetry of the initial and the scattering 
data. 

Statement 6.1. Initial data q{x,0) = q(x,0,e) given by ([7]) is even iff' 



a(-x) = -a(x) . (86) 
for all x in the domain of analyticity T> ofaix). Here a (x) = a(x)+ib(x) = — ^S'(x)+iA(x). 
Proof. If q(x, 0) is even then a(x) is odd and b(x) is even. Then a(—x) = — a{x) + ib{x) = 



a(x) — ib(x)] = —a(x). Conversly, taking a; 6 I, we obtain a(— x) + ib{— x) = a(—x) 



—a(x) = —a(x) +ib(x). □ 
Statement 6.2. Equation ( IBS"]) is equivalent to 



x{-z) = -x(z) (87) 

that holds for all x G B. 



Proof, x = x(z) z = a(x) <t=> — z = —a(x). Then, according to ( 1861) . — z 
that x(—z) = —x(z). Proof of the converse statement is similar. □ 

Before proving the next statement, it's worth reminding that the radical R(x, z) with a 
fixed x has a branchcut 7(2;) = •y m connecting a(x) and a(x) = a{x) that passes through 
//+. 

Statement 6.3. If a(x) is an odd function or if a(x) satisfies fl86|) then 



R 2 (-x,-z) = R 2 (x,z) , 
for all Qz > and all x ET>. The converse is also correct. 
Proof. Equation flHBl can be written as 



[{z - a(x))(z - a(x))} = [(z + a(-x))(z + a(-x))} . (89) 
The proof follows from equating coefficients of linear and free terms (in z). □ 

Remark 6.4. It is easy to verify using the moment conditions for a(x) from [21] that odd 
a(x) would imply w(0) = and 

W ' iC) -d( = . (90) 



ICI 



Equation (IBB]) obviously implies 



R(-x,-z) = ±R(x,z) . (91) 

Based on our choice of the branchcut of R, it is clear that in the case of real x, z the sign in 
(19T]) is positive if \z\ < ji + and negative if \z\ > In the case x 6 R the branchcut •y(x) is 
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symmetrical with respect to R and intesect R only at z = If x G R and ^z > 0, the sign 
in (1911) is positive only if there exists a curve cr(z) connecting z with the origin, such that 
a(z) fl j(x) = and simultaneously — a(z) R j(— x) = 0. It is clear that if such a(z) does 
not exist then a(z) can be choosen in such a way that a(z) fl 7(2) = and —cr(z) intersects 
j(—x) only one time. Then the sign in (19~T]) is negative. In the case of a complex x G £> the 
sign in (19~T]) is defined as above. 

Let E^x) denote the region (not necessarily simple) bounded by the curve j(—x) U 
{— 7(x)}, and let £(x) = (E*(x) U {— E*(x)} flC + . Clearly, S(x) is symmetrical with respect 
to imaginary axis and has the segment [— //+,//+] a its lower boundary. Then, for every 
x E B, equation (IBB]) has positive sign if z G S(x) and negative sign if 2; E(x). 



Theorem 6.5. If a(—x) = —a(x) then 



h(-x,-z) = Th(x,z) , (92) 

where the sign is negative if z G E(x) and positive otherwise. Conversly, ( l9~2"j) implies that 
a(x) is odd or a(—x) = —a(x). 

Proof. According to ( |50|) and to Statements 16.21 16.31 

x(z) r x ( z ) r x ( z ) 



X 



h(—x,—z)= / R(u,—z)du = — / R(—y,—z)dy — ^ / R(y, z)dy = ^fh(x, z 



(93) 

where y = —it. The last to terms of (193]) have sign minus if z G S(a;) and sign plus if 
z £ £(x). 

To prove the converse, we need to show that ( 1931) implies ( 1881) . Let us choose some 
x = iy, y (zM., z E S(x) and let us introduce 

m(y, z) = —ih(iy, z) = —ih(x, z) . (94) 

Then 

h(—x, —z) = h{iy, —z) = im(y, —z) — h(x, z) = im(y, z) , (95) 
so, according to fl93|) . 



m(y, —z) = m(y } z) . (96) 

Since m is analytic in y, we have Taylor expansion m(y,z) = YlT=o a k(. z )y k 1 so that (19~6"1) is 
equivalent to 



a k (z) = a k (-z) \/k G N . (97) 

Since m y (y,z) = J2T=i k> a k(z)y k ~ l has Taylor coefficients satisfying (!97|) . we obtain 



But m y (y,z) = h x (iy,z) = h x (x,z) = —R(x,z). Thus, (198"]) implies (191]) with positive sign. 
By analyticity of h, this result can be extended from purely imaginary to all x. Similarly, 
we can obtain (!9~T]) with negative sign when z G" E(x). Thus, ( 192]) implies (188]) . Statement 
16.31 completes the proof. □ 
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Corollary 6.6. In the casex,z G R equation (|92|) implies thatw(z) = sign(z — /j, + )Qh(x, z) 
is an even function and 

x, — z) = — z) if \z\ < ji + 
$lh(-x,-z)=m{x 1 z) if \z\> n + . (99) 

Corollary 16.61 states that ^sh(x, z) is even on z G [— and is odd outside [— 
Therefore, h z (x,z) is odd on z G [— and is even outside [— Taking into 
account Remark 16.41 and Statement 16.11 we obtain the following corollary. 

Corollary 6.7. If^sh(0,0) ^ i/ien (192]) implies that the initial potential q(x,0) is an even 
function. 

Given a scattering data (see [21]), we know w(z) but not necessarily h(x,z). So, does 
even w(z) implies even q(x, 0) under the inverse scattering procedure of [22], |24j ? According 
to ( |50l) . /i(x, z) depends on the branchcut of R(x, z). Let us denote by Rr, Rl and Kr, the 
radical R(x, z) (and corresponding to it h(x,z)) with the branchcut passing through 
respectively (in the case of even w(z) we have //_ = —//+)■ In all the statements above, we 
considered i? = Rr and /i = h R . In the case of x, z G K we choose Rl(x, z) to be positive 
for 2 > — fi + and preserve the orientation of contour 7 + (from /i + to — in the matrix 
RHP for the inverse scattering transform. With such choice of Rl, it follows from (1501) that 
hj_,{x,z) = ±hji(x,z) for x,z G M and |z| < //+ or |^| > /i + respectively. If a(x) satisfies 
f l86|) . then, according to Theorem 16. 5[ h^—x^—z) = hR(—x,—z) = —h R (x,z) for z G S(a?) 
and hi{—x, —z) = — Hr{— x, —z) = —hR(x, z) for 2; G" S(x). For 16I that means 

Sh L (-x, -z) = Sh R {x, z) (100) 

for all z. Thus, we can formulate a stronger version of Corollary 16.71 

Corollary 6.8. If w(z), z el, is even and if for all x > the inverse scattering procedure 
ofJEBi, W$ produces a semiclassical solution (j7J) with a(x) = —^S'(x) andb(x) = A{x), such 
that a(x) is analytic in a region containg x > and satisfies fl86l) in a vicinity of x = 0, 
then q(x,0), evenly continued to the whole M, is the initial potential corresponding to the 
scattering data w(z) according to the procedure of [22], \2Jfl . 

Proof. Let a(x), x > 0, be solution of the modulation equations (moment conditions), 
expressed in terms of w(z), see [21], Sect. 3.1.1. By Remark 3.1 of [21], a(x) = —a(—x) 
satisfies modulation equations when x < 0. Thus, we can assume that a(x), x G K, is 
analytic on M. and satisfies (!86l) . According to Theorem 16.5] the corresponding Hr and 
hi satisfy (|92l) . By the assumption, the main and complementary arcs have the required 
distribution of signs of S/i^(x, z) for all x > 0. Then, according to (llOOp . in the case of x < 
the main and complementary arcs have the required distribution of signs of Qh^x, z). Thus, 
the initial data represented by a(x), x G R, indeed corresponds to the given scattering data 
w(z) through the inverse scattering procedure of [22], [21]- □ 
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7 Examples 

7.1 Example of [22], ji = 2 

The case of 



a(a;) = ^tanha;, b(x) = sech (x), (101) 
/i > 0, was studied in details in [22]. The leading order (in e) 

is 

fo(z) = lim — \nr init (z,e) , (102) 

e— >0 Z 

where ri n a(z, e) is the reflection coefficient of Zakharov-Shabat problem for the focusing NLS 
with the inital data (GO), A(x) = b(x) and — |S"(x) = a(x), is calculated there to be 



fo(z) =(| - z) 



Z7T U 

h n(- 

2 v 2 



+ ln(z + T) + - — - ln(z - T) 



7T 

Ttanh -1 — H — In 2 H — e, when 3^ > 0, 
f 2 2 ' ~ ' 



(103) 



where T = y ^ — 1. This calculation involved exact solution of Zakharov-Shabat problem 
in terms of hypergeometric functions, followed by Stirling's asymptotic formula. Note that 
fo(z) has a log point T G £ in the case \i < 2 and has no log points in £ in the case 
fi > 2. This is the reflection of the fact that the initial data f llOip is purely radiative in the 
case /i > 2, and contains points of the descrete spectrum (solitons) on the vertical segment 
[— T, T] in the case \x < 2 (see [20]). From (11031) one readily obtain 

fo{z) = -y - ln(| - s) + 1 ln(, 2 - T 2 ) , (104) 

and 

= --signz (l - X[-t,t](z)) , z G R. (105) 



Since the asymptotic inverse scattering transform, developed in [22], [23], uses only the values 
w'(z), zeR, and, in the case fi < 2, the jump of fb(z) over the slit [0, T], we will focus on 
calculating these quantitees. 

Consider first the simplest case \i = 2. Then fi± = ±1, 

sinhx + i sinhx — i 

a { x ) = 1 j a{x) = (106) 

cosh x cosh x 

and 

n ^ 1 — isinhx 

a'(z) = p . (107) 

cosh x 

According to (11061) . a(x) is a meromorphic function with poles at y + 27um and zeroes 



at — y + 2irim, m G Z. Since 



, . , Z7i\ cosh£ ± 1 ,„ n . 

"(^y^^r- (108) 
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i^(+oo) 



-oo(-l) 



0(i) 



(0) 



Figure 12: The strip B = {x : — | < < |} is mapped onto the upper half plane of 
2-plane by z = a(x), so that a(B~) = £ + . Shown in parentheses are the corresponding 
images of the map a(x), i.e., the preimages of the inverse to a(x) map x(z). 



we see that B is the strip — ~ < Qx < |, where y = x(±oo), — y = x(0) and ±00 = x(±l) 
respectively, see FigO Since a'(x) 7^ inside B, there are no branchpoints in B. Note that 
for all x G R we have 



a 2 (x) + b 2 (x) = 1 



(109) 



so that 8 is the unit disc. 

In the case fi = 2, the inverse function x(z) to a(x) on C + is very simple. In order to 
find it, we calculate 



[z - a)(z - a) = (z - a) 2 + b 2 = z 2 - 2az + 1 = 



(110) 



which yields solution a(z) = Mz + -}. (Here we use notation a(z) for a(z) = a(x(z)).) Thus 



x(z) = tanh 1 a(z) = tanh 1 - \z H — ] 



According to (17B . we have 



px /-tanh x 

sign(l— z)$s / \/ z 2 — 2z tanh y + ldy = sign(l— z)Q / \J z 2 — 2zr] + 1- 

JxO) ia(z) 1 



111) 



dr\ 



w[z) ...... 

— 7] z 

(112) 

where z G K, x G K is arbitrary and y = tanh 1 77. Let us consider the case z > 0. 
Then a(,z) > 1, so the contour of integration in the latter integral is along the real segment 
[tanhx, a(z)} of the complex 77-plane (complex a-plane), except the singular point 77 = 1, 
which is traversed from above or from below if z > 1 or z < 1 respectively, see Fig. [T3J Since 
the integrand is real on M, only integration around rj = 1 contribute to the imaginary part 
of the integral. Thus 



w[z) 



iri Res 



^z 2 -2277 + I 



V 



r,=l 



— 1 — z) 

2 v ' 



(113) 



where the square root is choosen positive when z < 1 and negative when z > 1 respectively. 
By similar argument (but without sign(l + z) involved), we obtain 



w(z) 



VK 
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1 + z) 



(114) 



1 



fl(z) 



tanh x 



Figure 13: Upper contour of integration is for z > 1, lower - for z < 1. 



for z < 0. Thus, w(z) obtained through (ITTj) with initial data from (IllOp coincides with 
^sfo(z), zg8, where /o is given by (11031) with /i = 2 (and T = 0). The correct distribution 
of signs of $sh(x, z) along the main and complementary arcs (inequality fl2~Tl) ) follows directly 
from the analysis of zero level curves of Qh(x, z) with a fixed x G M in the upper z-halfplane 
(three zero level curves entering C + at ±1 and oo respectively meet at the branchpoint 
a(x)). One can also observe that fo(z) given by (I103P can be obtained through the Cauchy 
transform of w(z) (see [23]), and that sign distribution (J27j) for this fo(z) was proven in [22J. 



7.2 Example of [22], general case 



We consider now the general case \i > in (jlOip . Then (11061) - (1 11 01) become 



77 sinh x + i # sinh x — i 

a(x) = , a(x) = . 

cosh x cosh x 



Tj-isinhx «tt fcosh^il 

a(g)= _._ 1 .2„ » «U±y) = 

4 



cosh 2 x 2 7 sinh ^ 

\a 2 {x) +b 2 {x) = 1, 
(z — a) (z — a) = (z — a) 2 + b 2 =z 2 — /iz tanh x + 1 + T 2 tanh 2 x = 

respectively. The latter equation implies 



(115) 



tanha;(2;) = 5 , (116) 

1 1 

where the correct branch of the radical in the right hand side should be chosen. Equations 
(11151) show that £ is the ellipse centered at z = with horizontal and vertical semiaxes ^ 
and 1 respectively. The complement of £ + in C + is mapped onto the strip < Qx < ? in 
the x-plane. 

To find the image of £ + , we need to study the branchpoints x* of x(z) defined by 

sinh x* = — ~i . (117) 

The images of branchpoints in z-plane (log points), obtained by substituting (I117P into the 
expression for a(x), are 

z* = ±\— -1 = ±T. (118) 
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2 v_ 




Figure 14: Mapping a(x), case \x < 2. Shown in parentheses are the corresponding images 
of the map a(x), i.e., the preimages of the inverse to a(x) map x(z). Set B is the strip 
{x : — | < Qx < |} without the region bounded by the dashed contour, which is the image 
of two sides of the slit [0, T] by the map x(z). 




+ w(\ill2) 



(-T) i% (T) 
2 



Figure 15: Mapping a(x), case > 2. Shown in parentheses are the corresponding images 
of the map a(x), i.e., the preimages of the inverse to a(x) map x(z). Set B is the strip 
{x : — | < Qx < |} without the region bounded by the dashed contour, which is the image 
of the segment [— T, T] by the map x(z). 



In the case /i > 2 log points ±T 6 (— | , |) are not in the upper halfplane. In the case 
fi < 2 log points ±T = ±«|T|. In the latter case, we make a branchcut of x(z) in C + over the 
segment [0, T]. It is now easy to find the image B~ of the region S + (with the cut if \i < 2), 
see Fig. [TH Fig. [15] for /i < 2, /i > 2 respectively. Equations of the arc, connecting points 
x(0 ± ) and x(T), case /i < 2, and x(±T) and x(0), case /i > 2, are obtained from conditions 
3?a(?/) = and = respectively. These equations are 

— sin 7] = — cosh £ and cosh£ = — sin 77 (119) 

respectively, where y = £ + ir). 

Let us now calculate w(z), z G M, for the case \i > 2 (we only consider |z| > T) or for 
the case \i < 2. Using (I115p . we obtain direct and inverse maps between the values of a and 
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z by 

z = a+ijl-±a', a=2i^EZ!.| (120 ) 

The latter maps S + into a part of the lower a-halfplane and C + \S into the upper a-halfplane, 
see Fig. HH Fig. TH Using ffT2"UD . it is easy to check that if z > T, fi > 2 or z > 0, \i < 2, 
then the corresponding a = a(z) = a(x(z)) e K and a > ^. Making change of variables 
a(j/) = 77 in 



w(z) = sign(^ - z)9 / v/(z-a(y)) 2 + 6 2 (y)d?/ (121) 
^ A( 2 ) 

and using (I115p . we obtain 



/i /" W / 4r/ 2 drj 

w(z) =sign(- - / i (z-r]) 2 1 1 




2 



Direct calculation show that the integrand is real-valued. Thus, a contribution to the imag- 
inary part of the latter integral can come only from integration around the singular point 
r] = 7j. Repeating the previous argument (case \i = 2), we obtain 

«>(*) = £(£=F*) (123) 

if signz = ±1 respectively. Using the above technique, one could check that in the case /1 > 2, 
w(z) should be constant on the segment [— T, T] C K. Thus, w(z) satisfies the inequalities 
(|26|) for all \i > 0. In the pure radiational case /i > 2, the proof of the sign distribution 
(I27|) is the same as for the case of 11 = 2 in Section 17.11 This proof can be extended to the 
case m < 2 providing that the contour 7 does not intersect the branchcut of fo(z). However 
Lemma 4.11 from [22] shows that 7 m = r y m (x), which is defined by fo(z) given by (I103p . 
intersects the branchcut [0,T] if x is negative with sufficiently large \x\. In this case one 
can use Corollary 16.81 from Section [6] to show that q (x,0,e) is 0(e) close to q(x,0,e) for all 
x e JR. 

wzmgenpm 

Let us now calculate the jump 

Af(z) = f(x, z_) - f(x, z+) = h(x, z+) - h(x, z-) (124) 

for the case < ji < 2, where the points z± are equal but located on the opposite shores of 
the oriented vertical segment [T, 0]. According to fl50|) . we have 

Af(z) = T (2+) ^(z-a(y)y + V(y)dy (125) 

Jx(z-) 

where x(z-), x(z + ) are symmetrically located points on the dashed curve, Fig. [141 Using 
the change of variables rj = tanha;, we obtain 

Af(z) = [ V+ V^-^ + l + TV ■ , (126) 

J v- IT] 
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(1/2 



I174T 



-iW2T 



Figure 16: Mappings (11201) in the case fi > 2. The dashed line in the a-plane is the ima 
of R by the map o(z) = a(x(z)). In particular, d(±oo) = ±oo, a(±^) = a(±T) 

i^p, d(0) = 




(i 2 /4T 



Figure 17: Mappings (11201) in the case /x < 2. The dashed line in the a-plane is the ima 
of K by the map a(z) = a(x(z)). In particular, d(±oo) = ±oo, d(±^) = d(0 ± ) 

2 

±^r, d(T) = The solid line shows the contour of integration from a_ to a+. 
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where, according to (I120p . 



V± 



%z ± Vz 2 - T 2 Uy ± J-y 2 - T 2 



J^2 



2^2 



(127) 



with z = iy, see Fig. [TD, where a± = ^r]±. 

Observing that r]± are zeroes of the radical in (I126p . we obtain 



A/(z) = i ! v /z 2 -iiz V + l+T 2 V 2 -- 

J A 



(128) 



where A C C is a closed, clockwise oriented curve that contains the segment [r]_ , r] + ] . Thus, 
Res \/ z 2 — [izri + 1 + T 2 r] 2 ^ ± Res \/ z 2 — fxzr] + 1 + T 2 r] 2 

m{z — T) . 



ri=±l 



VK 

~2 



z 2 -fxz + ^-Jz 2 + fiz + ^-2T 



(129) 



Thus, we obtained the same value of Af(z) as (11031) has. It is now possible to reconstruct 
fo(z), given by (flQ3l with fi < 2, from u;(z) defined by ([123) and A/(z) defined by (11291 on 
[0,T]. 



7.3 The "Y" - shape of Bronski 

In [2], J. Bronski studied numerically discrete spectrum of ZS problem ([T]) with the potential 
g(x, Oe) given by <£T§, where 

A(x) = sech 2x, S(x) = /xsech 2x. (130) 

He found that for real valued q(x, Oe), i.e., for fi = 0, the accumulation curve for discrete 
eigenvalues of ([1]), (I130j) in C + in the limit e — > is a segment [0, i] of the imaginary axis. 
In the case \i — 1, the accumulation curve has "Y" - shaped form, with 2 = located at the 
bottow of "Y" . The change of shape of the accumulation curve happens at the critical value 

3 

fi* = 2~2. In this section we show that the endpoints of the accumulation curve coicide with 
the branchpoints of fo(z), obtained from the potential ([7]), (I130p by the AH transformation 
(T51|). In particular, /i* is a critical value when one branchpoint of fo{z) splits into two. 
For potential ((Tj), (I130p . we have 

. . usinh2a; i a sinh 2x + i cosh 2x /^„,s 

a(x) = 5 1 = — . (131) 

cosh 2x cosh 2x cosh 2 2x 



Then 



cosh 2 2x — 2 sinh 2 2x sinh 2x a cosh Ax + i sinh Ax — 3u 

a\x) = 2(j, 5 2% 2 — = ~- 3 > ( 132 ) 

cosh 2x cosh 2x cosh 2x 



so, after some algebra, equation a'(x) = for the branchpoints can be written as 

(fi + i)e 8x - 6/ie 4x +(i-i = 0. (133) 
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Solution of (I133P is given by 



v ' l > 1 ji + 1 

Thus, for fi = we obtain e 4x = 1, so x = is a branchpoint and z = a(0) = i is the 
corresponding log point in the spectral plane. 

The critical value fi* is given by equation 8fi 2 — 1 = 0, which yields fi* = 2~z. For 
[/,<[!*, direct calculation show that \e 4x \ = 1 for solution e 4x given by (11341) with the 
positive sign. Thus, the corresponding log point x G iR. Substituting fi* into (I134p . we find 
e 4x * = |(1 — 2y/2i), where x* = — | tan _1 (2-\/2) is the the branchpoint that corresponds to 
fi*. Direct calculations show that the corresponding log point on the spectral plane (the 
double point on Fig. 11 in [2]) is given by z* = a(x*) = ~ 0.91856z. 

To calculate the log point (s) (endpoint(s) of the accumulation curve) for other values of 
/i, we rewrite H133[) as 

{fi 2 + 1) cosh 2 4x - 6/i 2 cosh4x + 9fi 2 - 1 = , (135) 

which yields 



(cosh4x) 12 = 271 ' (136) 



3fi 2 ± y/1 - 8fi 2 

fi 2 + l 

Substituting fi — 0, we see that for fi < fi* we need to choose the positive sign in 
Calculating now cosh 2a; and sinh2x, we can express logpoints a(x) for an arbitrary fi > 
as 



-fi J 2fi 2 -l±Jl- 8fi 2 + id Afi 2 + 1 ± J 1 - 8fi 2 I 

a(x) = — V V 1 V^TT) , (137) 

4/i 2 + 1 ± a/1 - 8fi 2 

were only the plus sign should be used for fi < fi* (note the choice of signs for the radicals 
representing cosh2x and sinh2a;). 

Expression (I137j) provides the endpoints of the two "legs" of the "Y" -shaped accumulation 
curve (as e — > 0) for the points of the discrete spectrum in the case fi > fi*. In the case 
< fi < fi*, it gives the tip of the vertical segment on the imaginary axis where the points 
of the discrete spectrum accumulate. 



7.4 Double hump initial data 

Another example of an interesting initial data is 

a(x) = tanhx, b(x) = sech x — fcsech 2 x, k G [0, 1], (138) 

which contains double hump cases. Indeed, if k > ~, then b = u — ku 2 , where u = sech x, has 
a global maximum at u = l/2k G (0, 1). So, x = ± cosh _1 (2/;;) are two points of maximum 
of b(x). 
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Figure 18: Intersections of functions 1 + sin?] and 2fctan7? and points 773 



The points of ramification of the map 

2 ,, . 1 — i sinhx + 2zfctanhx 

a[x) = tanhx + «(sech x — fcsech x) with a (x) = 5 (139) 

cosh x 

satisfy equation 

1 — i sinhx + 2ik tanhx = 0. (140) 

Coefficients of this equation are 27ri periodic functions. The substitution u = sinh x in (11401) 
yields a forth order polinomial equation in u. Thus, there are no more than 4 ramification 
points of a(x) within a horizontal strip of width 27r in the complex x-plane. 
Separating real and imaginary parts of (11401) . we obtain the system 

I sinh(2£) cos(2??) — 2k sinh £ cos rj — sinh £ sin rj = (141) 
I cosh(2£ ) sin(2??) — 2k cosh £ sin r\ + cosh £ cos r\ — 0, 

where x = £ + irj. The first equation in (11411) has a common factor sinh£. Setting it zero, 
we obtain £ = 0. Then the second equation (I14ip becomes 1 + sin?? = 2/ctan??. Fig. [TBI 
shows that this equation has one positive solution 773 and one negative solution r}\ on (— n, it), 
where < rj^ < | and 7/4 < — |. We denote by Xg 4 the corresponding ramification points in 
the complex x-plane. 

Considering £ 7^ 0, we rewrite the first equation in (11411) as 

cosh£ cos(2??) — 2k cost] — sin?/ = 0, (142) 

which yields 

sin?? + 2A;cos77 

cosh£ = — . 143 

cos(2?i) 

Substitution of (I143j) into the second equation of (I14ip yields 

1 (4k 2 - 1) sin(2??) - 2k cos(2??) + \ S m(2 V ) cos(2??) 

cosh £ = — — r-= . (144) 

2 sm(2r/) cos(2r/) 

From (I151l) - (ll44p after some algebra we obtain 

u 3 + Ak 2 u + Ak = 0, (145) 



34 



where u = sin(2^). Equation (11451) should yield the remaining two points of ramification. 
In the case k = \ equation (11451) becomes [u + l)[u 2 — u + 2] = 0, which has the only real 
root u = —1. Substituting the corresponding r] = — | into the second equation of f 1 1 4 1 j) . 

we obtain cosh 2^ = 2\/2 cosh£, or £ = ±cosh _1 ^1 + -j^j- Thus, for k = \ we obtain the 
remaining ramification points 

4 2 = ±cosh- 1 (l + ^-^. (146) 

It is clear that in the case of arbitrary k > equation (11451) has only one real solution. 
Indeed, if that is not the case, then (I145P would have a multiple real root for some k > 0. 
But equation m'(u) = 3u 2 + k 2 = 0, where m(u) = u 3 + Ak 2 u + 4k has no real roots, so 
there is a unique real root u(k) of ( 11451) . Moreover, u(k) G (0,-1], since m(0) > and 
m(— 1) < 0. Thus, 7] = |sin _1 n(A;) or r\ = — | — |sin^ 1 n(/c). Substituting one of these 
values into (11431) (the former if k > | and the latter if k < |), we find the real components 
of the two remaining points of ramification x\ 2 . Note that ^tx\ = — tflx^- 

To calculate the pre- images of z = of the map (11391) . we set equation 

z cosh 2 x = sinh x cosh x + i cosh x — ik . (147) 

For the same reasons as above, for every z G <C + this equation has four roots in the strip 
—7i < Qx < 7r. Assuming z G K. and separating real and imaginary parts of (I147p . we obtain 



z(cosh(2£) cos(2r/) + 1) = sinh(2£) cos(2r/) — 2 sinh £ sin 7] 
2;(sinh(2£) sin(2?7) = cosh(2£) sin(2?7) + 2 cosh £ cos 77 — 2k 



'148) 



If z — then the first equation yields: 1) sinh£ = 0, or; 2) 

cosh£ cos(2?7) = sin 77 (149) 
Substituting £ = in the second equation (11481) yields 

cos 77(1 + sin 77) = k, (150) 

which has exactly one positive and one negative root on (— | , | ). The latter root is denoted 
x(0), see Fig. M 

To analize the remaining roots of (11471) with z = 0, we substitute 

cosh£ = ^f (151) 

COS 277 

into the second equation (I148j) with z = 0. After some algebra, we obtain 

sin 3 2t? + 2k sin 2 2t? - 2k = . (152) 

Introducing u = sin2r7, we can rewrite (I152p as 

m{u) =u 3 + 2ku 2 - 2k = . (153) 
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Figure 19: The graph of 1±^2 
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Figure 20: The image C and the domain B of the map z = a(x), where a(x) is defined by 
(I139j) . Ramification points x*, i — 1,2,3, 4 in the x-plane and the corresponding log points 
in C are marked by "crosses". Vertical cuts in C and their pre-images in B are shown by 
dashed lines. Round brackets in the x-plane are used to denote pre-images of points. The 
region C consists of the upper z-halfplane with three cuts. The upper boundary of B consists 
of two curves, connecting the pre-images of ±1 with y (pre-images of z > 1 and z < — 1 
respectively), and the closed dashed curve (the pre-image of the cut from a(x|) to zoo). The 
lower boundary of B consists of two dashed curves that contain points x^x^ (pre-images 
of verical cuts in £ + ), the curve that connects them and passes through x(0), and curves 
connecting the dashed curves with ±oo (the latter three curves form a pre-image ob the 
interval (—1, 1)). 
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It is easy to see that m(u) has a real root u(k) G [0, 1). Let us show that m(u) has only 
one real root. Indeed, m'{u) = (3u + 4k)u. So, critical points of m(u) are Ui = — 4^ and 
m 2 = 0. It is clear that U\ is a local maximum, w 2 is a local minimum and m(-u 2 ) < 0. Since 
m{u\) = (p/c 2 — 2) k < for all k G [0, 1], we conclude that m(u) has a unique real root 
u(k) and u(k) G [0, 1). 

Thus, we obtain sin27/ G [0,1), which means 2t] G [0, 7r] or 2r/[— 2ir, — tt]. That means 
rj G [0, |] or r] G [— 7r, — f ]. However, according to (11511) . we have a restriction 



which imply that 



sin 77 

7T > 1, 154 

cos 277 

sin 7] < — or — < sin i] < —= . (155) 

v 2 2 y 2 

Combining (I155p with the above restrictions on 77, we obtain 

3-7T 7T 7T 7T , , 

< 77 < or — < 77 < — . (lob) 

4 ' 4 6 - ' 4 v ; 

However, we will show that latter option is not possible. Indeed, implicitly differentiating 
(11531) . we obtain 

du 2fl - u 2 ) 

> 



dk (3m + Ak)u 

for any u G (0,1). Thus, the root u{k) is monotonically increasing. Substituting u(k) = 
sin 277 = sin | into (I153p . we obtain the corresponding k = > 1. Thus, for all k G [0, 1], 
the corresponding u(k) < ^, so that (11541) cannot be satisfies for the corresponding 77. 

Thus, we proved that for any k G [0,1] there exists a unique 7/ and unique cosh£, 
connected through (I15ip . that satisfy ()148p . Moreover, such i] G (— f)- Thus, we 
proved existence of the second pair of roots of a(x) = located in the strip —^f< ^sx < — f 
and symmetrical with respect to the imaginary axis. 

Obtained information allows us to sketch the domain B for the map ( I138p . see Fig. [20j 
Calculation of w(z) and of the jumps Afo(z) over the branchcuts is not included into this 
paper. 
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